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We introduce a solid on solid lattice model for growth with conditional 
evaporation. A measure of finite size effects is obtained by observing the time 
invariance of distribution of local height fluctuations. The model parameters 
are chosen so that the change in the distribution in time is minimum. On 
a one dimensional substrate the results obtained from the model for the 
roughness exponent a from three different methods are same as predicted 
for the Kardar-Parisi-Zhang (KPZ) equation. One of the unique feature of 
the model is that the a as obtained from the structure factor S(k, t) for the 
one dimensional substrate growth exactly matches with the predicted value 
of 0.5 within statistical errors. The model can be defined in any dimensions. 
We have obtained results for this model on a 2 and 3 dimensional substrates. 
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The KPZ equation [1] is one of the most studied stochastic equation in 
the field of growth. 

^ = u V 2 h + X(Vh) 2 + r ] (1) 
at 

Here, h(r,t) is the height function, A is coupling parameter, and r](r,t) is 
Gaussian noise with the correlation < r)(r,t)r)(r',t') >= 5(r — r',t — t'). 
First term is a linear term [2] referred as Edward- Wilkinson (EW) term. 
Applications and uses of KPZ equation have been well demonstrated [3]. In 
particular its use in understanding growth phenomenon has lead to vigorous 
activities in development of theoretical methods [4, 5], lattice models [6, 7, 8] 
and numerical methods [9] built around Eq. (1). The critical exponents 
in (1+1) dimensions are exactly obtainable [10]. However, in all the higher 
dimensions the determination of these exponents has been a difficult task. It 
is well known that this equation shows phase transition in dimensions higher 
than its critical dimension (2+1) [10] as a function of its coupling parameter. 
For weak coupling case the coupling constant renormalizes to zero leading to 
a linear equation. For strong coupling case, the perturbation approach fails 
and other methods are required. However, obtaining exact values for the 
exponents has never been possible although the ranges in which exact values 
are expected to fall are evident from the available references [4, 5, 6, 7, 8, 9]. 

Various lattice models have been devised [6, 7, 8] that are known to be- 
long to the KPZ universality in the asymptotic region. Most of the models 
suffer from finite size effects arising from cutoff length a and the substrate 
size L. One of the tests to probe the presence of finite size effects is to deter- 
mine the growth exponents by different methods such as height-height (h-h) 
correlations, structure factor S(k,t) =< h(k,t)h(—k,t) >, from saturated 
widths W sa ti etc.. For any given lattice model, the values obtainable from 
these methods can be statistically different. Since all these models are ex- 
pected to converge asymptotically to KPZ behavior, the apparent mismatch 
of the exponent values from different methods will be due to the finite size 
effects. On the other hand if a model gives same exponents within statisti- 
cal errors, it is expected to be free of finite size effects. In the following we 
propose a model that we believe to belong to the KPZ universality and, in 
(1+1) dimensions, it provides the values of the exponent a same within the 
statistical error, using three different methods of determination. This value 
also compares well with the exactly known value in (1+1) dimensions, 0.5. 

We describe the model below and the changes therein for (2+1) and (3+1) 
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dimensions. A site is chosen randomly and height at the site is increased by 
unity signaling random deposition on the substrate. The deposited atom is 
conditionally accomodated, otherwise evaporated. In (1+1) dimensions the 
deposited atom is accommodated if both its neighbors have at least same 
height as the deposited one. Otherwise, largest of the height differences at 
the site % and the nearest neighboring sites ,Sd = max(hi — hj),j = i + — 
is obtained and accommodation is allowed according to the probability factor 
e' s d/( 2a ). Thus Sd is the largest local step. Choice of a depends upon the 
behavior of the model for the given value of a. We choose the value of a 
that leads to minimum variation in the local height fluctuations. Model with 
such a a is expected to be least affected by the finite size effects [11]. It 
has been shown in reference [11] that a measure of finite size effects for a 
given lattice model can be obtained from the distribution of local height 
fluctuations. In this method, we define the local height (hi(t))i oca i = hi(t) — 
+ h i+ i(t))/2 with respect to the local reference as the average of 
nearest neighbor heights. Similarly we measure (hi(t + At))i oca i where At > 
w(t). w(t) is the width of the interface at t, and the inequality ensures that 
the difference between local heights measured at t and t+ At are uncorrelated. 
Thus we measure the distribution of uncorrelated fluctuations Ah(t)i oca i from 
the difference Ah(t)i ocal = {hi(t))i ocal - (hi(t + At))i oca i . Fig. 1 shows such a 
distribution for o = 1.7. Another distribution is obtained at later time and 
compared with the earlier one. In our case we have obtained distributions 
at * =500 MLs and t =5000 MLs for comparison with At = 100 MLs. Since 
the counts at Ah(t)i oca i = are largest in the distribution, the statistical 
error is minimum for zero fluctuation. We therefore use the parameter P — 
100( " f500 i ~^ 5000 ), where It is the count at Ah(t)i oca i = to measure the time 
invariance of the distribution of Ah(t)i oca i in (1+1) and (2+1) dimensions. 
In (3+1) dimensions, At and times for comparison are smaller due to the 
large computation times involved. Ideally P should be zero. In the present 
context we look for a minimum value of P as a function of a. The ratios 
Pq are obtained by averaging over large enough runs so that values of Pq 
are statistically discriminated for different values of a. It has been shown in 
reference [11] that this method is useful in identifying presence of finite size 
effects for any lattice model belonging to KPZ or EW universality. Fig. 2 
shows the variation of Pq with a for the model described. We have measured 
Po for the model with different a values on a substrate of length L = 40000. 
As can be seen, the minimum occurs at a = 1.7. We have therefore used this 
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value in (1+1) dimensions. For other values of a we found that value of a as 
obtained from structure factor deviates from 0.5 and the linear range is also 
reduced on the log- log plot. This confirms the effectiveness of the method 
for determining finite size effects [11]. 

In (2+1) dimensions the deposited atom is accommodated if three or more 
nearest neighbors have at least same height as its own. If this condition is 
relaxed to less number of in-plane neighbors, cross over due to EW region 
is obtained. The cross over is negligible when direct accommodation with 
three or four in-plane neighbors is allowed. For depositions at the site with 
less than three in-plane neighbors, the accommodation is decided from the 
largest of the four steps around the site using the exponential probability 
factor e~ s d/( 2 °" ). In (2+1) dimensions we have observed that a = 2.5 shows 
minimum P = 0.002 ± 0.0015%. The substrate size is L = 400. 

In (3+1) dimensions the deposited atom is accommodated if 5 or 6 nearest 
neighbors have at least same height as the deposited atom. We have used a = 
4.5 in the simulations since this value gives minimum Po — 0.0012 ± 0.001%. 
The At = 20MLs and the distributions are compared for the times 50 MLs 
and 500 MLs for the substrate size of L = 100. 

We present results obtained from (h-h) correlations, W sat as a function of 
L , and structure factor. The (h-h) correlation is 

x' 

where, correlation length £(t) ~ t l l z . In the limit x — > 0, / — ► 1. 
the width over a substrate of length L as, 

w 2 (L,t) = ^(h(x,t)-h(t)) 2 

X 

Here h(t) is the average height at time t. It can be shown that [10] for times 
t » L, g — > const. , thus w sat oc L 2a . For t « L w 2 (t) oc t 2(S . 
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The structure factor S(k, t) is measured as 

S(k,t) =< h(k,t)h(-k,t) > (4) 

where h(k,t) = (1/JV) £ x (/i(x,£) ~ MM))e ik - x . 

Figures (3), (4) and (5) summarize the results obtained for the model in 
(1+1) dimensions. Fig. 3 shows the log-log plot of G(x, t) Vs. x. The straight 
line is fitted between x =8 and 500. The slope gives a = 0.504 ± 0.002. Fig. 
4 shows plot of W sat Vs. L for (1 + 1) dimensions. The least square fit to the 
points gives a = 0.502 ± 0.005. 

Fig. 5 shows log-log plot of < h(k,t)h(—k,t) > Vs. k. The straight line 
fit is between /c=0.03 to 1.25. The slope near k = ir tends to zero [12]. The 
slope is 2.003+0.021. This gives [10] a = 0.500 + 0.021. Earlier, in reference 
[13] for the etching model, a slope 1.92 ±0.02 was obtained in the range of 
k=0.05 to 0.1, resulting in a =0.46. This was considered to be one of the 
best value obtained for the existing lattice models by this method. Clearly 
the present model provides a better value. In the same reference, a = 0.496 is 
obtained from the W sat . The apparent difference in the two a values indicates 
the presence of finite size effect for the etching model. In the present model, 
the slope is unaffected at smaller k values. We have tested it up to k = 0.005. 

Above results show that the proposed model has minimal finite size ef- 
fects. It further confirms the method introduced in reference [11] for determi- 
nation of finite size effects in a lattice model. We have applied this method in 
(2+1) dimensions and in (3+1) dimensions to choose the model parameters 
corresponding to the minimal finite size effects. 

In Fig. 6, we display the results of W sat Vs. L in (2+1) dimensions. 
The filled squares are the values calculated from the simulation results. The 
corresponding fit gives a = 0.357 ± 0.005. Fig. 7 shows the log-log plot of 
G(x,t) Vs. x. From its slope, a = 0.355 ± 0.001. The line is fitted between 
x = 2 to 50. Thus both the methods give results matching within statistical 
margin. 

In Fig. 6, we have also plotted the results of W sat Vs. L for the model in 
(3+1) dimensions. It gives a = 0.289 ± 0.005. The a values obtained from 
these models in (2+1) and (3+1) dimensions are very close to those predicted 
in reference [14]. In this reference a = 0.35702 for (2+1) dimensions and 
0.28125 for (3+1) dimensions. 

We have also measured the /3 values for these models from the log-log 
plots of w 2 (t) Vs. t. We obtain (3 = 0.332 ± 0.001 ,0.221 ± 0.002 and 
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0.168 ± 0.003 for (1+1), (2+1) and (3+1) dimensions respectively . These 
values are consistent with the universal relation a + z = 2 for the KPZ 
equation [10]. 

In conclusion, we have developed a lattice model belonging to the KPZ 
universality with minimum finite size effects. That the finite size effects are 
minimum is evident from the results obtained. The a value using W sat Vs. 
L , G(x, t) Vs. x , and from S(k, t) Vs. k plots are equal within statistical 
margin in (1+1) dimensions. The simulation values are very close to the 
exact a, 0.5. Measurement involving local height fluctuations is successfully 
used in determining the finite size effects in lattice models . The models 
with minimum finite size effects are expected to lead to the better accuracy 
in determining the exact exponents for KPZ growth in higher dimensions. 
In (2+1) dimensions, we have obtained values close to a — 0.36 while in 
(3+1) dimensions it is around 0.29. Both these values are close to the earlier 
prediction in reference [14]. 
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Figure 1: Plot of distribution of Ah(t)i oca i for the (1+1) dimensional model 
discribed in the text on semi log scale. The distribution is for t = 5000 MLs, 
with a = 1.7. The distribution is obtained by collecting the data over 3000 
runs. 



Figure 2: Plot of P in % as a function of parameter a for the model in (1 + 1) 
dimensions. 



Figure 3: Plot of G(x,t) Vs. x on log-log scale. The substrate size L=80000, 
and the number of monolayers grown are 5xl0 5 . 



Figure 4: Plot of W sat Vs. L on log-log scale for the model in (1+1) dimen- 
sions. 



Figure 5: Plot of S(k,t) Vs. k on log-log scale for the model in (1+1) 
dimensions. The points are averaged over the substrate lengths of L=800, 
900,1000,1100,1200,1300,1350. k varies from tt/100 to ir. 



Figure 6: Plot of W sat Vs. L on log-log scale for the model in (2+1) dimen- 
sions (filled squares) and for the model in (3+1) dimensions (filled circles). 



Figure 7: Plot of G(x,t) Vs. x on log-log scale for the model in (2+1) 
dimensions. The substrate size is L = 800, and the number of monolayers 
grown are 20000. 
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